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I. INTRODUCTION 



The inclusion of gravity into a physical theory fundamentally alters its basic assumptions 
and structure. Perhaps the most general and universal of all theories is thermodynamics. 
It is therefore particularly interesting to understand the changes induced by gravity in its 
underlying principles and to organize them into a general framework that serves as the 
foundation of a thermodynamic theory that incorporates strongly self-gravitating systems. 
Although this is relevant in its own right, it is also important in a somewhat different 
context: gravitational thermodynamics is expected to arise as one of the macroscopic limits 
of quantum gravity. Despite the fact that progress continues to be made into the way gravity 
alters the structure of quantum field theories [||], there does not exist yet a complete theory 
of quantum gravity. We believe that the characteristic features and principles of a theory 
of gravitational thermodynamics have to be understood properly before they can be fully 
justified by means of statistical methods derived from one or several candidate theories of 
quantum gravity. 

It is our purpose in this paper to formulate the general principles of gravitational thermo- 
dynamics. We shall discuss in detail the basic definitions and concepts of this formalism, its 
overall structure, its fundamental problem, the minimal set of assumptions (the postulates) 
that lead to the formal resolution of this problem as well as their mathematical and physical 
consequences. This formulation clarifies not only the differences and similarities between 
gravitational and nongravitational thermodynamics within a coherent framework, but also 
a number of existing misconceptions concerning its logical foundations and results. It gen- 
eralizes the overall structure of thermodynamic theory and may provide a basic framework 
to incorporate current and future progress in a statistical mechanical description of 
self-gravitating systems. 

A clear and elegant formulation of the foundations and structure of ordinary (that is, 
nongravitational) thermodynamics has been introduced by Callen 0. However, this formal- 
ism cannot describe the thermodynamics of gravitational systems in its present form and 
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needs to be modified. To recognize this, we shall follow the following strategy. We believe 
that the most effective way to gain insight into the limitations of a set of physical principles 
is by studying a model system that illustrates them without unnecessary complications. 
We introduce such a model problem in Sec. II. It is the simplest example of a composite 
self-gravitating system at finite temperature and resembles as closely as possible the text- 
book example of a composite system in ordinary thermodynamics. The model is general 
enough to capture all the non-trivial behavior of gravity, but simple enough to allow an exact 
evaluation of all quantities, and provides insight into the places where the thermodynamic 
formalism has to be refined. As a preparatory step in the study of the postulatory basis 
of thermodynamics, we show that it provides a counterexample to a basic assertion of 
ordinary thermodynamics: in gravitational thermodynamics, additivity of entropies applied 
to spatially separate subsystems does not depend on or require the entropies of the latter to 
be homogeneous first-order functions of their extensive parameters. 

The general principles of gravitational thermodynamics are developed from the start in 
Sec. Ill, which constitutes the core of the paper. It is shown that a general and rigor- 
ous characterization of the defining properties of thermodynamical self-gravitating systems 
originates from two factors, namely (1) the particular characteristics of their extensive vari- 
ables, and (2) the homogeneous properties of their intensive parameters as functions of 
extensive variables. The correct definitions of these variables are discussed in detail. Ex- 
tensive variables include quasilocally defined quantities like quasilocal energy and provide 
the background for the postulates. We discuss the composition of thermodynamic systems, 
the existence of fundamental equations, and formulate the fundamental problem of gravita- 
tional thermodynamics. Our object is to follow as closely as possible the logic of Callen's 
formalism and generalize its basic definitions and postulates wherever it is required. The 
definitions and postulates arrived at are a natural extension of the ones of ordinary ther- 
modynamics with appropriate modifications necessary to accommodate the global aspects 
and nonlinearities characteristic of gravity. The postulates are the basic principles of ther- 
modynamics when strongly self-gravitating systems are considered. We revise completely 
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the logic of Sec. II and show that additivity of entropies (and, in general, additivity of all 
Massieu functions) as well as the generalized first and second laws of thermodynamics are 
part of these postulates and, as such, are not to be proven within the thermodynamic formal- 
ism. We demonstrate that the fundamental problem is formally solved by these postulates 
in spite of strong interaction between constituent subsystems and the associated nonaddi- 
tivity of extensive variables. The preeminence of Massieu functions over thermodynamic 
potentials is discussed. The conditions of thermodynamical equilibrium and its (spatially) 
inhomogeneous nature (equivalence principle) are a direct result of the postulates. 

Although the fundamental equation and the intensive variables maintain their mutual 
relationships, the former is no longer a homogeneous first-order function of extensive param- 
eters. This property is not a consequence of the inhomogeneity of equilibrium configurations 
but of the functional form of the gravitational equations of state. These are in general no 
longer homogeneous zeroth-order functions under rescaling of extensive parameters. We 
show that these central differences with nongravitational thermodynamics are not forbidden 
by the logical structure of the formalism and can be incorporated easily into it by relaxing 
several assumptions in the postulates. However, the mathematical consequences of the new 
postulates are quite different from the ones familiar in regular thermodynamics. Formal 
relationships (like the Euler equation) must be reformulated and there is no direct analogue 
of the ordinary Gibbs-Duhem relation among intensive variables. Nonetheless, this does not 
influence signiflcantly the applicability of the thermodynamic formalism. 

The principles of the framework and its logic apply to general systems at finite temper- 
ature. They incorporate the so-called gravothermal catastrophe and other peculiar thermo- 
dynamical behavior observed in self-gravitating objects. As such, all standard equilibrium 
thermodynamics of gravitational configurations is a consequence of the general postulates 
presented here. The approach generalizes ordinary thermodynamics and provides a consis- 
tent treatment of composite self-gravitating systems (with or without matter fields) . More- 
over, it becomes evident that the modifications in the thermodynamic formalism necessary 
to incorporate gravity liberates it from assumptions that appeared (and are not) fundamen- 
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tal, and as such, allows us to regard its general character and the extent of its logic in their 
full measure. 



We motivate our presentation of the postulates of gravitational thermodynamics and 
their consequences with a model problem. It consists of a spherically symmetric uncharged 
black hole surrounded by matter (represented by a thin shell). We shall find in this section 
the relationship among the entropy functions 5*0, Sb, and Ss whenever thermodynamic 
equilibrium is achieved and its connection with the functional dependence of those functions. 
In what follows, subscripts B and 5* refer to quantities for black hole and shell constituents 
respectively, whereas the subscript refers to quantities for the total system. 

The composite system is characterized by the surface area Aq = AttTq of a two-dimensional 
spherical boundary surface Bq (located at r = tq) that encloses the components, and the 
quasilocal energy Eq contained within ^,0. The Arnowitt-Deser-Misner (ADM) mass m+ 
of the system as a function of these variables is 



Throughout the paper, units are chosen so that G = fi = c = kBoitzmann = 1- Due to 
spherical symmetry of the model problem, we use areas and radii interchangeably. The 
pressure po associated to the surface Bq is obtained as (minus) the partial derivative of the 
energy Eq with respect to Aq at constant m+. It corresponds to a negative surface tension 
and reads 



where kQ = (l — 2m+/ro)^^^. Let Pq denote the inverse temperature at the surface Bq. The 
first law of thermodynamics for the system reads 



II. HOMOGENEITY AND ADDITIVITY 




(2.1) 




(2.2) 



dSo = (3q {dEo + Pq rf^o) , 



(2.3) 
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which can be written as a total differential by using Eqs. ( |2.1|) and ( p.2|) as 

dSo = (3o ko~^ dm+ . (2.4) 

Consider now the constituent systems. The black hole is characterized by the surface 
area Ab = AttR"^ of a spherical boundary surface Bji located at r = i? < tq, and by its 
quasilocal energy Eb- The ADM mass m_ of the black hole as a function of these variables 

isi 

m4EB,AB)=EB(^l-^^ . (2.5) 

The horizon radius of the black hole is 2m__. The different radii satisfy the inequalities 
2m„ < 2m+ < -R < Tq, where 2m_|_ represents the horizon radius of the total system. The 
pressure associated to the gravitational field of the black hole at the surface Bji is 

MEs.As) ^ ^(l - f ^ ^(1 - t.f , (2.6) 

where fc_ = (1 — 2m_/i?)^/^. As it is well known, the first law of thermodynamics for the 



black hole can be expressed as a total differential by using Eqs. ( |2.5| ) and (|2.6| ), namely 

dSB = Pb [dEB + Pb dAB) 

= i3Bk_-^dm_, (2.7) 

where (3b is the inverse temperature of the black hole at the surface Bji. This equation must 



be contrasted with the entropy differential (|2.4| ) for the total system. 



For thermodynamical purposes, the shell is considered (effectively) at rest. For simplicity, 
we assume that the shell surface coincides with the surface Bb- (This does not imply loss 
of generality 0].) Thus, the surface areas for the black hole and shell coincide: Ab = As = 
An = AttR'^. The gravitational junction conditions at the shell position determine its 
surface energy density a and surface pressure ps to be 

Es = AixR^a = R{k.- k+) , (2.8) 
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and 

PS = \k-(l - k+f - k+il - A;_)2 

respectively. The symbol Es denotes the local mass-energy of matter and A;+ = (1 — 
2m+/i?)^/^. To simplify the analysis further, we consider only the case when the total 
number of particles A^^ in the shell is constant. The condition > m_ guarantees that 
both Es and "ps are non- negative. Let denote the inverse local temperature at the shell. 
Its entropy differential reads 

dSs = i5s{dEs+VsdAn) ■ (2.10) 

Use of Eqs. ( p.8| ) and ( |2.9| ) allows us to write the matter entropy as 

dSs = Ps {k+~^ dm+ - k_-^ dm_) . (2.11) 

How are the three entropies related when the system is in equilibrium? The emergence 
of equilibrium conditions from general principles within thermodynamics is the subject of 
the following section. However, intuitively the system is in a state of thermal equilibrium 
provided (1) (3b = Ps = Pr, and (2) (Sr = (3q ko~^ fc+ at the surface Br. The first condition 
constrains the temperature at the shell surface to coincide with the black hole temperature 
there (black hole and shell in mutual thermal equilibrium), whereas the second guarantees 
that the total system is in thermal equilibrium with its components 0. Mechanical equi- 
librium of the matter shell with the black hole is guaranteed by the shell pressure equation 
Under these conditions, Eqs. (U), {^), and (^TTO|) jointly imply that 

dSo = dSB + dSs . (2.12) 

The entropy of the composite system is, therefore, additive with respect to its constituent 
subsystems. Since the entropy is a function of energy and size (its "extensive variables" 
discussed below) this means that, up to a global constant, 

So{Eo,Ao) = SB{EB,An) + Ss{Es,An) . (2.13) 
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(2.9) 



In the preceding analysis, additivity is a direct consequence of the conditions of thermal 
and mechanical equilibrium. Observe that it is independent of the functional form of the 
parameters Ps{Es, Ar) and [3b{Eb, Ar). This is as expected, since inverse temperature 
appears in the first law of thermodynamics as an integrating factor. In particular, the fore- 
going derivation of additivity does not depend on the special choice of boundary conditions 
or phenomenological matter action employed in Ref. ^ or on spherical symmetry [l^. It 



depends only on the adopted definition of stress-energy tensor JlTl in terms of quasilocal 
quantities. 

As discussed in the following section, the entropies Sq{Eq,Aq), Sb{Eb, Ar), and 
Ss{Es, Afi) can be determined from Eqs. ( |2.4| ), (|2.7| ) and (|2.10|) only if the precise forms of 



all the thermodynamical equations of state are known. The latter express intensive param- 
eters as functions of the appropriate extensive parameters. For example, it is well known 
that if the thermal equation of state for a black hole is given by Hawking's semiclassical 



expression W2 



MEb, An) = SttEb (l - (l " ^) ^ (2-14) 



then Eq. ( |2.7|) yields the Bekenstein- Hawking entropy 

SBiEB,AR)=A7TEB^(^l-^'j =A7imJ. (2.15) 

These expressions are sufficient to demonstrate that additivity of entropies for spatially 
separate subsystems does not require the entropy of each constituent system to be a ho- 
mogeneous first-order function of its extensive parameters. This is in contrast to ordinary 
thermodynamics 0,^. Recall that a function f{xi, a;„) is said to be a homogeneous m-th 
order function of the variables (xi, ...,x„) if it satisfies the identity 

/(Axi,...,Ax„) = A"^/(xi,...,x„) , (2.16) 

where A is a constant. Upon the rescaling Eb ^Eb, Ar — > X^Ar {R XR) the entropy 
Sb{Eb, Ar) in Eq. ( ^.15] ) behaves as a homogeneous second-order function of Er and as a 
first-order function of Ar, namely 



Sb{XEb, X'An) = \' Sb{Eb. Ar) . (2.17) 



Equations ( |2.14] ) and ( p.6|) also illustrate a central characteristic of gravitational systems: 



the inverse temperature and pressure are not homogeneous zeroth-order functions. Un- 
der rescaling they behave as [3b{XEb,\^Ab) = X Pb{Eb, A^) and pb{XEb, X^Ab) = 
X-^Pb{Eb,Ar). 

The functional form of Ss{Es, Ar) and its homogeneous properties depend on the explicit 
form of the matter equations of state. These arise from either a phenomenological or a field 
theoretical description of the matter fields involved, and their precise form does not concern 
gravity. Examples of equations of state for a self-gravitating matter shell (in the absence of 
a black hole) in thermal equilibrium with itself have been studied in Ref. [|l3l. Observe that 
Eq. (|2.2| ) for pq and Eq. (|2.6| ) for pb are indeed equations of state while Eq. (^^) for ps is 
not. If the equations of state were known for both components, the total entropy Sq could 
be obtained by Eq. ( p.l3|) as a function of the extensive parameters of the subsystems. A 



discussion of this point and of further properties of this model are delayed to the following 
section. 



III. THE FUNDAMENTAL PROBLEM AND THE POSTULATES 

The preceding analysis motivates the search for principles that are independent of model 
problems and that incorporate the characteristics of gravity into a logical framework more 
general than ordinary thermodynamics. We must start, therefore, by reviewing the basic 
assumptions. 

Gravitational thermodynamics describes (effectively) static states of macroscopic finite- 
size self-gravitating systems. As in regular thermodynamics, it is expected that very few 
variables survive the statistical average involved in a macroscopic measurement. What are 
these macroscopic parameters? The relationship between thermodynamical and dynamical 
variables in gravity has been amply discussed and we refer the reader to the literature |T^-|T^. 
For our purposes, it is enough to recall the following points. Firstly, it has been shown in 
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a wide variety of problems (involving black holes at finite temperature in interaction with 
matter fields) that the thermodynamical energy coincides with the quasilocal energy E that 
naturally follows from the action of a spatially bounded region. If denotes the three- 
dimensional boundary of the system and the two-surface resulting from its intersection 
with a spacelike hypersurface S, the quasilocal energy is the value of the Hamiltonian that 
generates unit time translations on ^B in the direction orthogonal to the surface ^B |ir^ , [TT| . 
We postulate that this is the appropriate energy variable in (gravitational) thermodynamics 
for all self-gravitating systems. Secondly, the analog of the thermodynamical "size" of the 



system is the induced two-metric aab of the two-dimensional boundary surface B [y,|T8|,|l4 
This property reflects the "surface character" of gravitational thermodynamics and is in 
part a consequence of the lack of an operational definition of "volume" in the presence of 
black holes. The size reduces to the surface area A of the two-surface only in the case of 



spherical symmetry |TJ]. For composite systems, quantities that measure size for internal 
matter components have to be found. Finally, the remaining macroscopic variables are a 
finite number of conserved charges. These may include, for example, angular momentum 



TM, suitable combinations of electric [113 or magnetic pO] charges, cosmological constant 



other types of hair |^|, and number of particles for matter systems |]T3[. (The ther- 



modynamic conjugate quantities to these parameters are chemical potentials defined at the 
boundary of the system ||T9| , p!^ .) 



The existence of these macroscopic parameters motivates the first postulate: 
Postulate I: There exist particular states (called equilibrium states) of self-gravitating 
systems that are completely characterized macroscopically by the specification of a finite set 
of variables. These variables are the quasilocal energy, size, and a small number of conserved 
quantities (denoted generically by the symbol A^). 

In ordinary thermodynamics a similar postulate is usually applied to so-called "simple" 
systems . These systems do not include gravitational or electromagnetic fields and are by 
definition macroscopically homogeneous. The previous postulate incorporates strongly self- 
gravitating configurations (with or without matter fields). As shown below, these systems 
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may be spatially inliomogeneous. Systems describable by these parameters may be termed 
"simple" in gravitational thermodynamics. 

Observe that the preceding postulate does not imply that every gravitational system has 
equilibrium configurations. Very often a system does not possess an equilibrium state even 
though it has definite values of energy and other parameters. Rather, the postulate maintains 
that equilibrium states, in case they exist, are completely described by the foregoing finite 
set of parameters. 

The variables {E, aab, N) that describe a gravitational equilibrium state are to be called 
extensive parameters. Extensive quantities are constructed entirely from the dynamical 
phase space variables. Another essential difference with usual thermodynamics appears 
here: in the latter, the extensive parameters of a composite system equal the sum of their 
values in each of the subsystems. As we illustrate below, this is not the case in gravitational 
thermodynamics. 

Some extensive variables of a self-gravitating system cannot be constrained in the con- 
ventional thermodynamic sense. For example, there exist no walls restrictive with respect 
to angular momentum of a stationary black hole system. However, this is not unusual or 
particular to gravity. For instance, it also occurs in the treatment of magnetic systems: 
there exist no walls restrictive with respect to magnetic moment. However, one can main- 
tain always the value of magnetic moment constant at a boundary surface that delimits the 
system by a feedback mechanism that continually adjusts the magnetic field [Q. The same 
happens in gravitational thermodynamics: unconstrainable quantities can be kept constant 
at a given boundary surface by continually monitoring the value of their respective conjugate 
chemical potentials at this surface. The unavailability of walls that restrict certain extensive 
variables is only an idiosyncrasy that does not affect the applicability of thermodynamics. 
As in ordinary thermodynamics, a boundary that does (does not) allow the flux of heat 
can be called diathermal (adiabatic). Observe that the quasilocal energy adopted here has 
a very important property for thermodynamics. It is "macroscopically controllable" in the 
usual thermodynamic sense: it can be "trapped" by restrictive boundaries and "manipu- 
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lated" by diathermal ones. A boundary that does not allow the flow of heat and work can 
be called "restrictive with respect to quasilocal energy." A closed system is deflned as one 
whose extensive variables (quasilocal energy, size, etc.) remain effectively constant at its 
boundary surface. 

It is of course difficult to split a self-gravitating system into independent "component" 
systems in the manner familiar in ordinary thermodynamics. Although one can speak of a 
"composite" system formed by two or more subsystems, the latter interact strongly among 
themselves. Clearly, if a composite system is closed, the simple systems are not necessar- 
ily so. However, internal constraints may exist among the component systems. These are 
constraints that prevent the flow of energy or any other extensive parameter among subsys- 
tems. For example, in our model problem internal constraints can restrict the flow of energy 
between the two subsystems (for instance, by flxing a particular value of Eb) or area (by 
fixing Afj). The relaxation of internal constraints in an equilibrium composite system will 
create processes that will tend to bring the system to a new equilibrium state. 

The central problem of thermodynamics of strongly self-gravitating systems is, therefore, 
a reflection of the central problem of ordinary thermodynamics namely: The determi- 
nation of the equilibrium states that will result when internal constraints are removed in a 
closed, composite system. 

What assumptions are needed in order to solve this problem? Equilibrium states in 
gravitational thermodynamics must be characterized by a simple extremum principle. As any 
other thermodynamic system, a gravitational system will select, in the absence of constraints, 
any one of a number of states, each of which can also be realized in the presence of a suitable 
constraint. Each of these constrained equilibrium states corresponds to particular values of 
the extensive parameters of each constituent system and has a deflnite entropy. Therefore, 
the extremum postulate states that if constraints are lifted, the system will select the state 
with the largest entropy. Paraphrasing Callen: 
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Postulate II: There exists a single-valued function (the entropy 5") of the extensive vari- 
ables of any composite system, defined for all equilibrium states, and possessing the following 
property. In the absence of internal constraints, the values assumed by the extensive pa- 
rameters are those that maximize the entropy over the manifold of constrained equilibrium 
states. 

This postulate has to be interpreted carefully. Physical equilibrium states correspond to 
states that extremize the total entropy over the manifold of constrained equilibrium states. 
Equilibrium states are, therefore, either maxima, minima, or inflection points of the entropy. 
However, in the absence of constraints, the extensive parameters of the components in the 
final equilibrium state will be those that maximize the entropy. Postulates I and II not 
only predict equilibrium states, but also determine their stability properties. Equilibrium 
states corresponding to maxima of entropy are stable whereas unstable equilibrium states 
correspond to extrema other than maxima. It is important to emphasize that Postulate 
II makes no reference to nonequilibrium states. Furthermore, it implies neither that all 
equilibrium states of a gravitational system must have maximum entropy nor that stable 
states do exist. After all, it is common to find systems possessing equilibrium states that 
are local minima of entropy. Simple examples include a nonrelativistic self-gravitating gas 



in a spherical box or isothermal stellar systems |23| . 

The entropy as a function of its extensive variables constitutes the "fundamental equa- 
tion" of a self-gravitating system. The first differentials of the fundamental equation define 
its intensive variables. For systems with a vanishing shift vector, the intensive variables in 
the entropy representation are (/?, Pp, where P denotes inverse temperature, p pressure, 
and fi chemical potentials. Systems possessing a nonvanishing shift require functional dif- 
ferentials in the definitions of their intensive parameters. (This happens, for example, in the 
thermodynamic description of stationary geometries fl^-) For general spacetimes the con- 
jugate quantities to the size aab are proportional to (minus) the spatial stresses introduced 



in Ref. []TT|. The intensive parameters are always functions of the extensive parameters. The 



set of functional relationships expressing intensive in terms of extensive parameters are the 
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thermodynamical equations of state of a self-gravitating system. For example, for static (as 
opposed to stationary) systems these are 

(5 = f5{E,A,N) , 
(3p = (3p{E,A,N) , 

Pfi = Pfi{E,A,N) . (3.1) 

As in ordinary thermodynamics, once the fundamental equation S{E, aab, N) of a system 
(or, alternatively, its complete set of equations of state) is known, all its thermodynamical 
information can be obtained from it. 

The criteria for global and local stability of equilibrium states in terms of the entropy 
function are identical to the ones familiar in ordinary thermodynamics PJl3| . In particular, 
global stability requires that the entropy hypersurface S{E, aab, N) lies everywhere below 
its family of tangent hyperplanes. 

It is possible to express the fundamental equation in terms of different sets of inde- 
pendent variables by performing Legendre transformations on the entropy. These are the 
so-called Massieu functions |2^J§]- They play a more fundamental role in gravitational 
than in ordinary thermodynamics because they are in a one-to-one correspondence with 
actions [^]. Their preeminence over "thermodynamic potentials" has not been empha- 
sized sufficiently. (The latter are Legendre transforms of energy and include the Helmholtz 
potential F and Gibbs potential G.) For static systems, Massieu functions include, for ex- 
ample, Si{P, A,N) = S — 13 E = —PF, in which quasilocal energy has been replaced by its 
conjugate entropic intensive parameter (inverse temperature) as independent variable, and 
(S2(/3, Pp, N) = S — PE — PpA = —PG in which, in addition, the size of the system is replaced 
by its entropic intensive parameter. The above equations can be generalized to stationary 
geometries if one recalls that in general it is not possible to choose all intensive variables 



constant at a given choice of two-dimensional boundary surface |p!^ , p5| . The basic extremum 
postulate is very general and can be reformulated in these alternative representations: each 
Legendre transform of the entropy is a maximum for constant values of the transformed 
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(intensive) parameters 0]. 

It is important to emphasize that the second postulate incorporates not only the so-called 
first law, but also the generalized second law [p^ -|28[] into a thermodynamic formalism. 

Finally, how is the entropy of a composite system related to the entropies of the sub- 
systems? In the previous section we illustrated that entropies are additive despite strong 
interaction between subsystems. On the one hand, additivity of entropies does not seem to 
depend critically on the particular functional form of the intensive parameters. On the other 
hand, it seems to be a natural consequence of the additivity of actions in a path integral 
approach to statistical thermodynamics |]5|Jl0|]. These reasons motivate us to assume the 



additivity postulate: 

Postulate III: The entropy of a composite system is additive over the constituent subsys- 
tems. Furthermore, it is a continuous, different iable, and monotonically increasing function 
of the quasilocal energy E. 

We emphasize three important points. First, we shall not assume in this postulate 
that the entropy of each subsystem is a homogeneous first-order function of the extensive 
parameters. The postulate is, therefore, more general than the corresponding one in regular 
thermodynamics Second, the preceding postulate implies that all Massieu functions are 
additive over component Massieu functions. As we illustrate below, this is not the case for 
thermodynamic potentials |Q. Third, the monotonic property implies that the temperature 
is postulated to be non-negative as in ordinary thermodynamics. 

The logic of the previous section must be contrasted with the present one: additivity is 
neither the result of equilibrium conditions among intensive variables, nor of the functional 
form of intensive or extensive parameters, but a fundamental assumption. Additivity is valid 
even when the subsystems cannot be considered independent and strongly interact among 
themselves. As we show in the following paragraphs, equilibrium conditions are indeed a 
consequence of the postulates. 

The preceding postulates are the natural extension of the postulates of nongravitational 
thermodynamics necessary to accommodate the extensive parameters characteristic of grav- 
is 



itational systems. Are these postulates sufficient to solve the fundamental problem despite 
strong interactions among subsystems? The answer is on the affirmative. To illustrate this 
consider again our model problem in the light of the logic resulting from the postulates. We 
shall determine the equilibrium state of the closed, composite system, namely, the relation- 
ships that must exist among extensive variables of the subsystems for the total system to 
be in thermal and mechanical (and in general, chemical) equilibrium. We also shall indicate 
how far one can proceed in the explicit solution of this problem without assuming particular 
expressions for the equations of state of the subsystems. 

By Postulate I, the black hole and matter subsystems are simple systems characterized 
by the extensive variables {Eb,Ab,Nb) and {Es,As,Ns), respectively. The composite 
system is itself a simple system and is characterized by the variables (£^o; ^o, -^o)- The 
size of all systems reduces to the area of their respective surfaces. By Postulate II, the 
fundamental thermodynamical equations in the entropy representation are the functions 
So = So{Eo, ^0, A^o), Sb = Sb{Eb, Ab, Nb), and Ss = Ss{Es, Ag, Ng). Postulate III states 
that Sq{Eq, Aq, Nq) = Sb{Eb, Ab, Nb) + Ss{Es, As, Ng). For simplicity and with no loss of 
generality, we assume Ab — As = Ar and the quantities A^O; ^b and Ns to be constant. 
The system is considered closed if its quasilocal energy and area are kept effectively constant 
at the boundary Bq, namely 

Eq = const.; Aq = const. (3.2) 

The fundamental problem is to determine the extensive variables {Eb, Eg, Afj) as functions 
of these constant quantities whenever equilibrium is attained as a result of relaxing internal 
constraints. Postulate II establishes that the total entropy of a composite system in a state 
of equilibrium is an extremum, namely, it does not change as a result of an infinitesimal 
virtual transfer of heat or work from one subsystem to the other. Therefore, in equilibrium 

o?<S'o = = dSB + dSs 

= (3b {dEB + Pb dAn) + [dEs + ps dAn) , (3.3) 
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where the second equahty is a consequence of Postulates I and II. The entropic intensive 
variables are defined in the conventional way: 

(3s{Es,An) = {dSs/dEs)A^, PsPs{Es,Ar) = {dSs/dAR)Es; MEb,Ar) = {dSB/dEB)A^, 
and I3bPb{Eb,Ar) = {dSB/dAR)EB. 

Since the quasilocal energy can be expressed as Eq = r^il — ko), the closure equations 
are equivalent to the condition m+ = const. Because the energy Eb refers to the surface Br 
which coincides with the shell surface, it is easy to see that the total quasilocal energy at 
Br is 

Er = Eb + Es 

= R{l-k+). (3.4) 



This equation is a consequence of the additivity of quasilocal energy discussed in ||5|,|TT 
(If the black hole energy Eb is defined at a surface which does not coincide with the shell 
surface, the energies Eb and Es are not simply additive as in Eq. (|3.4| ) [^.) The closure 
conditions and Eq. ( p^) allow the total entropy (p.3|) to be written as 

dSo = = {/3s - /3b) dEs + {/3bPb + /3sPs - /3bPr) dAR , (3.5) 
where the pressure pr is defined as 

PH(^H.^«)-l|^(l-f)"'- 15^(1 -M^ (3.6) 



Since the equality in Eq. ( p.5[ ) must be satisfied by independent and arbitrary variations of 
Es and Ar, we must necessarily have 

(3s = /3b=(3r, (3.7) 

and 

Ps+Pb = Pr- (3.8) 

The preceding equations are the sought equilibrium conditions. They state the relation- 
ship among intensive variables of the subsystems for the composite system to be in thermal 

17 



and mechanical equilibrium. As in nongravitational thermodynamics, they yield a formal 
solution to the fundamental problem provided the equations of state 

/3b = I3b{Eb, Ar, Nb) , Pb = Pb{Eb, Ar, Nb) ; 

Ps = Ps{Es, Ar, Ns) , PS = PsiEs, Ar, Ns) (3.9) 

for the subsystems are known. If this is so, Eqs. (|3.7|) and ( p.8|) are two formal relationships 
among Eb, Es, Ar and m+ (with Ns and A''^ each held fixed). Equations (p^), ( p.4|) , (p.TD, 



and (|3.8|) are, therefore, the four desired equations that determine the four sought variables 
{EB,Es,AR,m+). 

Naturally, the variable m+ in Eqs. ( |3.2| ) and ( |3.4| ) does not play an important role in 
the formalism: the relationship among the three energies Eq, Eb, and Es can be written 
explicitly as 

For a closed system, Eqs. (|3.7|), (p^), and ( |3.10| ) (with Eq. (p.2|)) provide three desired 
equations to determine the three sought variables {Eb, Es, Ar). 

The fundamental problem is formally solved by the postulates despite the quasilocal 
energy Eq not being the simple sum of the component energies Eb and Es due to binding 
and self-energy interactions characteristic of gravitational systems. (It is easy to see, by using 
Eq. (|3.10|) that all thermodynamic potentials obtained from Eq by Legendre transforms are 
not the simple sum of the component potentials.) This indicates not only that the postulates 
form a complete set of assumptions for a more general class of thermodynamic systems than 
previously considered, but also the appropriateness of the adopted definitions of extensive 
parameters. 

Consider some further consequences of the postulates. Firstly, the mechanical equilib- 
rium condition ( |3.8| ) represents the spatial stress component of the junction conditions at 
the surface Br. (It reduces to Eq. ( p.9| ) if the pressure equation of state for the black hole 
is given by (^.6])). Secondly, additivity of entropies and the equilibrium conditions ( |3.7| ) and 
( p.8|) allow the differential of the total entropy to be written as 



dSo = Po {dEo + po dAo) 

= (3r {d{EB + Es) + {pB + Ps) dAR) 

= PnidEn+pRdAR) . (3.11) 

This expression confirms that there is no "gravitational" entropy associated to the shell 
29|,|5|| and illustrates that Er is the total quasilocal energy and pr the associated pressure 



associated to the surface Br. In turn, Eq. (|3.11|) implies that 

dSo = Po ko~^ dm+ = Pr k+~'^ dm+ . (3.12) 

Therefore, the inverse temperature at the surface Br is given in terms of the inverse tem- 
perature Po at the boundary Bq by 

Poh-' = PRk+-\ (3.13) 

The (spatially) inhomogeneous character of thermodynamic equilibrium (equivalence 
principle) is a consequence of the postulates of thermodynamics and the definition of quasilo- 
cal stress-energy. Thus, the postulates do incorporate equilibrium states in inhomogeneous 



systems in contrast to the ordinary postulates of thermodynamics where a system 
that is not homogeneous is not in thermodynamic equilibrium even if its properties remain 
constant in time. 

The preceding treatment of a composite self-gravitating system must be contrasted with 
the equivalent one of a composite nongravitational system presented in Appendix A. Al- 
though the systems are physically different, their similarities and differences are readily 
apparent. In particular, the gravitational equations ( p.lO| ) and ( |3.2| ) substitute the relations 
( |Al| ) and ([A2| ) of flat spacetime thermodynamics. 



The formalism provides the methodology to solve the fundamental problem for self- 
gravitating systems. In the spirit of thermodynamics, it yields explicit answers for explicit 
functional forms of the fundamental equations (or equivalently, the associated equations of 
state) of each of the subsystems These are outside the realm of thermodynamics and are 
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the result of either phenomenological or statistical mechanical descriptions of the constituent 
systems. We reiterate its formal structure: for a composite self-gravitating system one 
must assume the fundamental equation of the components to be known in principle. If 
the total system is in a constrained equilibrium state (characterized by particular values 
of the extensive parameters for each constituent system), the total entropy is obtained by 
adding the individual entropies of the components and is, therefore, a known function of 
their extensive parameters. The extrema of the total entropy define the equilibrium states. 
Stable states correspond to maxima of entropy. As an illustration, if explicit equations 
of state are known for both black hole and matter in our model problem, their entropy 
values can be found (up to an overall constant) by integrating Eqs. ( |2.7|) and (|2.1CI|) , and 



substituting the values of {Eb, Es, Aji) at equilibrium. The total entropy is then given by 
Eq. d^). 

Can this logical framework accommodate runaway instabilities (the so-called gravother- 
mal catastrophe) observed in bounded self-gravitating systems? The answer is on the af- 
firmative. This behavior is a consequence of the postulates and the particular form of the 
fundamental equations characteristic of gravity. Typically, the latter are such that there 
exist, besides equilibrium states that locally maximize the entropy, equilibrium states that 
locally minimize it. (The existence of these state is well-known in stellar dynamics [^] and 



black hole physics p,|T8|.) Consider as illustration an isothermal self-gravitating gas in a 
closed spherical container p3|. The system might be thought of as formed by two compo- 



nents, the 'core' and the 'halo.' The formalism states that if the entropies for the components 
are known as functions of their extensive parameters, the total entropy is S" = Sc + Sh- Equi- 
librium states are obtained by extremizing this function and are characterized by particular 
values of the extensive parameters of the components. For simplicity, consider only the ener- 
gies (or equivalently, the density contrast between components). The entropy functions for 
the gas components are such that there exists an equilibrium state (described by a particular 
critical value of the density contrast) that is a local entropy minimum over the set of all 
possible equihbrium states But Postulate II predicts that this state is unstable. The 
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onset of instability in the gas obeys the postulates: if the system finds itself in that state 
and the density contrast between core and halo is allowed to change, the system will try to 
reach equilibrium states of higher entropy. The system finds it advantageous to transfer en- 
ergy (or work) from one region to another, developing more internal inhomogeneities. Local 
stability conditions 0JT^ imply that a negative heat capacity is associated to a local entropy 
minimum: if the core gets hotter than the halo, heat flows from the core to the halo and 
the core temperature raises. The end result depends on the form of the entropy function 
and on the direction of the fluctuation that started the instability. It might be that a local 



entropy maximum exists in which the system settles down (as discussed in Ref. [|23[, this 
occurs if the entropy is such that, for example, Ch < \Cc\)- In this case the halo temperature 
rises more than the core's and the system shuts off in a stable state. A runaway instability 
happens if there does not exist a local maximum for the system to settle down (this happens 
if the fundamental equation is such that Ch > |Cc|)- In this case the temperature difference 
between halo and core keeps growing. Whether a black hole is created or the system runs 
out of equilibrium before that occurs, the important point is that, as long as the system can 
be described by equilibrium physics, the postulates predict its behavior if the fundamental 
equations of the components are known. The above argument applies equally to a collection 
of stars or other astronomical systems. 

We have studied so far the impact of gravitational extensive variables in the thermody- 
namic formalism. But the latter is also characterized by the functional forms of its intensive 
variables (equations of state). These arise from a dynamical theory but their main character- 
istic is that, in general, they are no longer homogeneous zeroth-order functions of extensive 
parameters (for instance, the inverse temperature Pb in Eq- ( |2.14]) is homogeneous first-order 
in energy and half-order in area; although the intensive variable jSp in the entropy repre- 
sentation of a static black hole is homogeneous zeroth-order, this is not the case for other 
systems [0.) This implies that the consequences of the postulates are different than in ordi- 
nary thermodynamics, particularly the mathematical properties of fundamental equations. 
Fundamental equations are in general no longer homogeneous first-order functions of their 
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extensive variables. (Alternatively, the homogeneous properties of fundamental equations 
in gravitational thermodynamics imply that intensive variables are no longer homogeneous 
zeroth-order functions.) This does not affect the formalism itself, but has direct implications 
for at least two formal relationships among thermodynamic quantities. Firstly, the so-called 
Euler relation is necessarily different from the one familiar in ordinary thermodynamics . 
An Euler relation is a consequence of Euler theorem stating that a homogeneous function 
/(xi, ...,x„) of m-th order satisfies the equality 

m/(.i,...,x„) = x,(^)+... + x.(|^) . (3.14) 

In standard thermodjTiamics entropy is a homogeneous first-order function and the Euler 
relation is therefore S = (3E + (3pV — (3^N. In contrast and as an example, the Euler relation 
for a static charged black hole reads 

S = ^(3E + (3pA-^(3fiN . (3.15) 

Euler relations for hollow self-gravitating thin shells with power law thermal equations of 



state have been presented in Ref. |13 



Secondly, there does not exist a Gibbs-Duhem relation in gravitational thermodynamics 
13| . The Gibbs-Duhem relation in ordinary thermodynamics is a direct consequence of the 



homogeneous first-order properties of the fundamental equation and relates the intensive 
parameters of a system. It states that the sum of products of extensive parameters and the 
differentials of the (conjugate) intensive parameters vanishes. In the entropy representation 
it reads Edf3 + Vd{f3p) — Nd{f3fi) = 0. In contrast, if one combines the first law with the 
Euler relation ( p.l5|) for a charged black hole one obtains 

Ed(3-(3dE + 2Ad{(3p)-Nd{(3iJ.) + (3^idN = . (3.16) 

The reformulation of an Euler relation and the absence of a Gibbs-Duhem relation set gravity 
apart from other interactions: even for magnetic systems the Euler and the Gibbs-Duhem 
expressions maintain their usual relationship. 
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There are no obstacles in applying the preceding formalism to any self-gravitating sys- 
tem. These may include not only nonrelativistic astronomical objects, but also highly rel- 
ativistic systems involving general black holes in interaction with matter. The (spherically 
symmetric) model problem was used only as an illustration because of its simplicity and 
transparency. For general situations, a larger state space is required to incorporate a larger 
number of extensive variables and thermodynamic equilibrium includes equilibrium under 
"interchange of number N." It is possible to employ a condensed notation where the symbols 
Xi and Pj denote generically all extensive and intensive parameters, respectively (exclud- 
ing energy and inverse temperature), as in Ref. 0. In this way the equilibrium conditions 
( p.7| ) and (|3.8|) are easily generalized to include all chemical potentials for the system. Al- 



though the treatment of these and other systems may be technically difficult, the resolution 
of the fundamental problem of thermodynamics obeys the same logical structure as the one 
presented here. 

We emphasize again that the characteristics of gravitational thermodynamics are the 
result of its extensive parameters (which must include quasilocal quantities like quasilocal 
energy) and the particular homogeneity properties of its intensive variables (equations of 
state) as functions of extensive parameters. This is a more general and rigorous character- 
ization of the defining properties of gravitational thermodynamics than, for example, the 



one presented in Ref. |^ in the context of black holes. The definitions of extensive and 
intensive variables as well as the changes introduced in Callen's postulates are the minimal 
changes necessary to incorporate the global aspects and nonlinearities characteristic of the 
gravitational interaction into a postulatory formulation of thermodynamics. 

To summarize, we have presented the overall structure and principles of a thermodynamic 
framework that incorporates self-gravitating systems. All the results of standard (gravita- 
tional) thermodynamics are a consequence of the generalized postulates and the solution of 
the fundamental problem and can be extracted from them by following the standard proce- 
dures described in Refs. P,|5^. Possible applications of the formalism include, for example. 



the description of quasi-static reversible and irreversible processes, alternative representa- 
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tions, phase transitions, etc. Finally, a further postulate is usually introduced in standard 
thermodynamics: the so-called third law. However, the main body of thermodynamics does 
not require this postulate since in the latter there is no meaning for the absolute value (and 
therefore for the zero) of entropy. The role and interpretation of this kind of assumption in 
the statistical mechanics of gravitation is the subject of a future publication PB| ]. 
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APPENDIX A: 

The simplest closed, composite system in ordinary thermodynamics consists of two sub- 
systems separated by a movable diathermal wall The fundamental problem is to 
determine the extensive variables {Ei,Vi, Ni) and {E2,V2, N2) of the subsystems whenever 
equilibrium is attained. The quantities E, V, and refer to internal energy, volume, and 
number of particles, respectively. The closure conditions are 

Et = Ei + E2 = const. , (Al) 
Vt = Vi + V2 = const. , (A2) 
Nt = Ni + N2 = const. (A3) 

The postulates of thermodynamics allow a formal solution of the problem. The conditions 
of equilibrium follow from the equation 

dST = = dSi + dS2 

24 



= Pi {dEi + pi dVi - fii dNi) + P2 {dE2 + P2 dV2 - fi2 dN2) 

= (A - (32) dEi + (A pi -(32P2) dVi - (A/xi - /?2/i2) dNi . (A4) 

The first equality is a consequence of tlie extremum value of entropy (Postulate II) whereas 
the second is a consequence of additivity of entropies (Postulate III). The last equality 
follows from the closure conditions. Since the variations dEi, dVi, and dNi are independent, 
Eq. (|A4|) implies that in the entropy representation 

(3i=(32, (A5) 
(3iPi=(32P2, (A6) 
f3i /ii =[32 1^2 . (A7) 

If the equations of state 

(3a = (3a{Ea,Va,Na) , 
Pa = PaiEa, K, A^a) , 

fia = fia{Ea, K, A^,) (A8) 

are known for the subsystems (with a = 1,2), then Eqs. ( |A5[ ), (|AB|), and ( |A7| ) provide 
three formal relationships among (£'i,Vi,A^i) and (£'2, V2, A^2)- Eqs. (|A1|) , (|A^), and (|A3|) 
together with Eqs. ([A5|), ([A6D and (|A7D are, therefore, the six desired equations necessary 
to determine the six unknown variables. 
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